Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN KTeX style file v2.2) 



A deeper X-ray study of the core of the Perseus galaxy cluster: the 
power of sound waves and the distribution of metals and cosmic rays 

J.S. Sanders* and A.C. Fabian 

Institute of Astronomy, Madingley Road, Cambridge. CBS OHA 



5 February 2008 



ABSTRACT 

We make a further study of the very deep Chandra observation of the X-ray brightest galaxy 
cluster, A 426 in Perseus. We examine the radial distribution of energy flux inferred by the 
quasi-concentric ripples in surface brightness, assuming they are due to sound waves, and 
show that it is a significant fraction of the energy lost by radiative cooling within the inner 75- 
100 kpc, where the cooling time is 4-5 Gyr, respectively. The wave flux decreases outward 
with radius, consistent with energy being dissipated. Some newly discovered large ripples 
beyond 100 kpc, and a possible intact bubble at 170 kpc radius, may indicate a larger level 
of activity by the nucleus a few 100 Myr ago. The distribution of metals in the intracluster 
gas peaks at a radius of about 40 kpc and is significantly clumpy on scales of 5 kpc. The 
temperature distribution of the soft X-ray filaments and the hard X-ray emission component 
found within the inner 50 kpc are analysed in detail. The pressure due to the nonthermal 
electrons, responsible for a spectral component interpreted as inverse Compton emission, is 
high within 40 kpc of the centre and boosts the power in sound waves there; it drops steeply 
beyond 40 kpc. We find no thermal emission from the radio bubbles; in order for any thermal 
gas to have a filling factor within the bubbles exceeding 50 per cent, the temperature of that 
gas has to exceed 50 keV. 

Key words: X-rays: galaxies — galaxies: clusters: individual: Perseus — intergalactic 
medium — cooling flows 



1 INTRODUCTION 

• The Perseus cluster, Abell 426, is the brightest galaxy cluster in 
the sky when viewed in t he X-ray band. Th e cluster contains a 

, bright radio source 3C 84 jPedlar et alJlT990h , the radio lobes of 
which are displacing the X-ray emitting thermal gas of the clus- 
ter (Bohring er et al. 1993; Fabian et al. 2000). The X-ray emis- 
sion from the intracluster medium (ICM) is highly peaked in the 
centre and the radiative cooling time of the hot gas is less than 
5 Gyr within a radi us of 100 kpc, decre asing to 3 x 10^ yr within 
the central 10 kpc jSanders et alj|2004h . A cooling flow of sev- 
eral 100 M0 yr^^ would take place if radiative energy losses from 
the inner ICM are not balanced by some form of energy injection. 
As expected from cooling, the gas temperature does drop from the 
outer value of 7 keV within the central 100 kpc but only down to 
about 2.5 keV with little X-ray emitting gas found at lower tem- 
peratures, except in coi ncidence with line-emittin g filaments seen 
at optical wavelengths jpabian et al.ll2003j . l2006l) . Heating by the 
central radio source is widely considered responsible for balancing 
the radiative cooling, although the exact mechanisms by which the 
energy is transported and dissipated and a heating/cooling balance 
established have been unclear. 
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Similar behaviour is found in ma ny X-ray peaked, cool-core 
clusters (see IPeterson & Fabiarj [20o3 for a review). The central 
galaxy in such clusters could grow considerably larger than ob- 
served if radiative cooling was unchecked. X-ray observations of 
their innermost regions provide an excellent means to study the 
feedback of an Active Galactic Nucleus (AGN) on its host galaxy 
in action. Here we use a long (900 ks) Chandra observations of 
the X-ray brightest cluster core to examine the energy balance and 
metallicity in considerable detail. The quality of the data, in terms 
of counts per arcsec, are much higher than available with any other 
cluster. We examine the power propagating through the core in 
terms of pressure ripples, or sound waves, and the cosmic-ray im- 
plications of a hard X-ray component. We also use the metallic- 
ity distribution, the temperature profile of the X-ray emitting gas 
across an optical filament and assess the thermal gas content of the 
radio bubbles. 

No X-ray emission has been detected from w ithin the radio 
bubbles jSchmidt et al] |2002| : ISanders et al] l2004h . Two depres- 
sions in the X-ray surface brightness are also found, to the north- 
west and south, not associated with high-frequency radio emis- 
sion. These are likely to be 'ghost' radio bubbles which have de- 
tache d from the nucleus and buoyantly risen in the gravitational 
field (IChurazov et ^l2000l : [pabian et alj|2000l) . an idea supported 
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by wea k low frequency radi o spurs seen pointing towards their di- 
rection jpabian et al.l2002bh . 

Surrounding t he inner radio lobes are X-ray bright rims 
jpabian et al. I l2000h at a higher pressure than the o uter regions, 
and s eparated from them by a weak isothermal shock jpabian et alj 
I2OO6I) . Extending further into the clusters are concentric fluctua- 
tions in surface brightness. These are plausibly pressure and density 
ripples, in which case the y are sound waves g enerated by the infla- 
tion of the radio bubbles ( iPabian et a l. 2003a, '2006"). The period of 
the waves is about 10 Myr, close to the expected age of the bubbles 
due to buoyancy, and is not plausibly related to any other source 
of disturbance. Rip ples due to s ound wav es have since been foun d 
in simulations (Rus zkowski et al. 2004; S iiacki & Springelll2006h . 
Such sound waves can transport significant energy in a roughly 
isotropic manner and so balance radiative cooling of the intra- 
cluster g as if they dissipate their energy over the surrounding 50- 
lOOkpc jpabian et alj|2003d; [Ruszkowski et aLll2004lFabian et al.l 
I2OO5I) . Concentric ripple-like features are also found around M87 
in the Virgo cluste r, and are interpreted there as weak shocks 
jporman et alj|200^ . More p owerful outbursts inclu de those dis- 
covered in MS 0735 .6-1-7421 jMcNamara et al.ll2005h. H ercules A 
jNulsen et al]|2005 j) and Hydra A i Nulsen et al]|2005bl) . 

Around the central galaxy in t he Perseus cluster, NGC 1275, 
lies a giant line-emitting nebula jLvndsl Il970l : IConselice et al.l 
l200lh. associated w ith cool ~ 0.7 keV X-ray emitting fila ments 
jpabian et al1l2003bh of ~ lO'' M© mass jpabian et al.ll2006h . The 



Ha emitting filaments appear to be have been drawn out of the 
central galaxy by the rising bubbles ( H atch et al..,2006. ). 

There is spectral evidence for hard X-ra y emission from the 
central region, f ound either with ho t thermal jSanders et alj|2004h 
or non-thermal jSanders et ai]|2005l) models. The 2-10 keV lumi- 
nosity of this emission is ~ 5 x lO^"' erg s^' . 

The iron metallicity structure in the core of the cluster is inho- 
moge neous and complex jSchmidt et alJl2002l : ISanders et alj2004l 
I2OO5I) with the abundance dropping in the very core. There are also 
structures such as high-metallicity ridge and blobs, which may be 
associated with the bubbles. 

We use a Perseus cluster redshift of 0.0183 here, which 
gives an angular scale of 0.37 kpc per arcsec, assuming Hq = 
70kms"' Mpc"'. 



ratio of its respective foreground dataset to the total foreground. 
The exposure time of each section was adjusted to give the same 
count rate in the 9-12 keV band as its respective foreground (where 
there is no source). The background sections were then reprojected 
to the original level 2 event file for the matching observation, and 
then reprojected to the 04952 observation. 

In addition we constructed separate background event files 
for each dataset to account for out-of-time events, where pho- 
tons hit the detector while it is being read out. We used the 
MAKE_READOUT_BG script (written by M. Markevitch) to con- 
struct these files from the original level 1 event files after filter- 
ing bad time periods. The readout backgrounds were reprojected to 
match the 04952 observation. 

Por the spectral analysis, for a particular region we extracted 
foreground spectra from each of the foreground event files rele- 
vant for the region in question. Similarly we made spectra from 
the background event files and and readout background files. The 
foreground spectra were added together to make a total foreground 
spectrum with a total exposure time. Background spectra were 
added similarly to make a total background, and so were readout 
background spectra. 

We also created responses and ancillary responses for each of 
the foreground datasets, weighting the responses according to the 
number of counts in each spatial region between 0.5 and 7 keV. 
A response and ancillary response was made for the total fore- 
ground spectrum by adding the responses and ancillary responses 
for the individual observations, weighting according to the number 
of counts between 0.5 and 7 keV. 

To analys e the spectra, we fit them in XSPEC version 11.3.2 
jAmauj[T99a) . The energy range 0.6 to 8 keV was used during 
fitting, and the spectra were grouped to have at least 20 counts per 
spectral bin. The spectra extracted from the appropriate region from 
the background files was used as a background spectrum, and the 
spectra from the region from the readout background files was used 
as a 'correction file'. 

In this paper we use the APEC |Smith et al.ll200lh and M EKAL 
jMewe et alJI 19851, 1 19861 : fKaastrall 19921 ; iLiedahl et alJI 19951) ther- 
mal spectral m odels. To model photoelectric absorption we use the 
PHABS model jBalucinska-Church & McCammonlll992l) . 



2 DATA PREPARATION 

The datasets analysed in this paper are those that were examined in 
|pabianetal.l ( [2006) . However here we used the standard CIAO data 
preparation tools on the event files. We used the datasets from the 
Chandra archive which had gone through Reprocessing III (repro- 
cessed using version 7.6.7.1 of the pipeline). Datasets 03209 and 
04289 were not reprocessed at the time of the analysis, so we there- 
fore went through each of the steps in the Science Threads to repro- 
cess the data manually for these. The reprocessing was done prior to 
CALDB 3.3.0 (using gainfile acisD2000-01-29gain_ctiN0005.fits). 

W e filtered the level 2 event files for flares as in Pabian et al. 
j2006h . We then reprojected each dataset to match the 04952 ob- 
servation in sky coordinates. We also reprocessed a combined 980- 
ks blank sky observation file to use the same calibration files as 
used by the foreground observations. We randomised the order of 
the events in the background file in order to remove any potential 
spectral variability. The background file was split into sections, to 
provide a background for each foreground observation. The length 
of each section was chosen to have the same ratio to the total as the 



3 X-RAY SURFACE BRIGHTNESS 
3.1 Surface brightness images 

IPabian et alj j2003j) and IPabian et alj j2006l) used unsharp- 
masking techniques to reveal the surface brightness fluctuations in 
the intracluster medium. Unsharp masking increases the noise in 
the outer parts of the image where there are relatively few counts. 
We have experimented with several techniques to improve on sim- 
ple unsharp masking. We split the exposure-map-corrected image 
(Pig.[T]left panel) into 40 sectors (centred on the central nucleus), 
and fitted a King model to each of the sectors outside of a ra- 
dius of 13 arcsec (to avoid the central source). We constructed a 
model surface brightness image by iterating over each pixel, us- 
ing the value obtained by interpolating in angle between the model 
surface brightness profiles of the two neighbouring sectors. This 
model image was then subtracted from the original image, resulting 
in Pig. [T] (centre panel). The image very clearly highlights the sur- 
face brightnes s increase associated with the low-tempe rature spiral 
in the cluster jpabian et alj200d ; IChurazov et alj|2000h . A number 
of other features can be seen, including the possible cold front to 
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Figure 1. Surface brightness images of the cluster. (Left) 0.3 to 7 keV full band X-ray exposure-map-corrected image, smoothed with a Gaussian of 1.5 arcsec. 
(Centre) Image after subtracting King model fits to 40 sectors, smoothed with a Gaussian of 1.75 arcsec. (Right) Original image after high-pass filtering, then 
smoothing with a Gaussian of 1 .5 arcsec. 




Figure 2. Full band X-ray image showing the sector examined for surface 
brightness fluctuations (Fig. [5). The region is centred on the inner NE radio 
bubble to better match the surface brightness contours. 

the south, the 'bay' and the 'arc' (figure 2 in iFabian et al.ll2006h . 
It does not show the ripples particularly clearly (except perhaps by 
the southern bubble) as the cool swirl is dominant. 

The ripples are more clearly highlighted using a Fourier high- 
pass filter technique. A two-dimensional fast Fourier transform of 
the exposure-map-corrected image was made. We removed the low 
frequency components with a wavelength greater than 75 arcsec. 
Frequency components between wavelengths of 75 and 38 arcsec 
were allowed through using with a linear filter increasing from to 
1 between these wavelengths. All shorter wavelengths were left to 
remain. The Fourier-transformed image was then transformed back 
to give Fig.[T](right panel) after light smoothing. This technique re- 
moves the cool swirl and the underlying cluster emission. It clearly 
reveals the ripples, presumably sou nd waves generated by the in- 
flation of the bubbles, discovered bv lFabian et alj i2003ah . It shows 
a previously unseen ripple near the edge of the image to the east. 

3.2 Sector across eastern ripples 

To quantitatively examine the ripples, we have examined the sur- 
face brightness in a sector shown in Fig.|2l We show in Fig.[3](top 
panel) a surface brightness profile (in the 0.3 to 7 keV band) made 
to the east-south-east (ESE) of the cluster core. Also shown in the 
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Figure 3. Surface brightness profile to the ESE of the cluster core (0.3 to 
7 keV), also showing a King model fit, fractional residuals from the fit, a 
temperature profile with fit, and fractional residuals from the temperature 
fit. The dotted lines on the temperature residual plot show the expected 
variation associated with a 3 per cent surface brightness fluctuation. Note 
that the profile is not centred on the nucleus. 
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top panel is a King model fit to tlie profile. In the second panel we 
show the fractional residuals to the fit, clearly showing the oscil- 
lations in surface brightness observed in Fig. [T] (right panel) at the 
few percent level (the first peak in this figure corresponds to the 
bright rim of the radio bubbles). 

To look for any temperature structure associated with the os- 
cillations we fitted PHABS absorbed APEC thermal spectral models 
to spectra extracted in annuli from the spectra. The fitting procedure 
allowed the temperature, metallicity, normalisation and absorption 
to vary, and minimised the C-statistic ( Cash 1 979) in each fit. The 
projected temperature profile is shown in the third panel of Fig. [3] 
We also show the results from fitting the same model (using the 
statistic) to deprojected spectra in larger bins. We used the depro- 
jection method in Appendix lAlto create the deprojected spectra. 

We finally s how the residuals from a simple 'rj model' 
jAllen et al]|200lh fits to the projected and deprojected tempera- 
ture profiles (Fig.[3]fourth panel) superficially. There is no obvious 
temperature structure associated with the ripples in this location 
of the cluster The amount of temperature variation expected from 
the surface brightness fluctuations can be estimated. Simulations 
of thermal spectra in XSPEC shows that the surface brightness is 
independent of temperature at constant density in the temperature 
range 3-6 keV. If the adiabatic index of the ICM is 7= 5/3, as- 
suming the ideal gas law and that the X-ray surface brightness is 
proportional to the density squared, the fractional temperature {T) 
fluctuations associated with surface brightness (/) changes should 
be of magnitude 



(1) 



8T _\8I 
~Y ~ 3T' 

We plot on the lower panel of Fig.|3]dotted lines showing the range 
of temperature variation expected to be associated with 3 per cent 
variations in surface brightness (which are the maximum observed 
here). These include a scaling factor of 5 to convert from projected 
surface brightness fluctuations to intrinsic emissivity fluctuations 
(see Section [33t . The deprojected temperatures are comparable to 
those expected but we caution that there is significant noise in the 
results. 



3.3 Wave power 

We now estimate the power implied by the surface brightness fluc- 
tuations observed in the cluster (Fig.[T]right), assuming that they are 
sound waves. This is then compared with the power radiated within 
the inner regions of the cluster core where the radiative cooling 
time is shorter than its expected age of a few Gyr (by age in this 
context we mean the time since the last major merger). 

The instanta neous power, transmi tted in a spherical sound 
wave IS given bv ( lLandau&LifshitJl959[) 



pc 



(2) 



at a radius r, where the sound wave pressure amplitude is 8P, the 
mass density of the medium is p and the sound speed is c. The 
sound wave pressure amplitude can be computed from the density 
amplitude with 



5 811 
SP= -n^kT—a, 
3 ne 



(3) 




assuming 7=5/3, where a is a factor to convert from the electron 
number density iig to the total number density. 

The fractional variation in density over the wave 8ng/ng is 



Figure 4. Power in surface brightness fluctuations, assuming they are due 
to spherical waves. The units are logjo erg s^' . The rectangles show the re- 
gions examined in Fig.|5]and|6] The circles show the excluded inner bubble 
and NW bubble regions. We assume a factor of 2.5 to convert from surface 
brightness to density variations. 



estimated from the surface brightness fractional change. As dis- 
cussed above, there is little variation of surface brightness with tem- 
perature at constant density in this temperature range. Assuming 
bremsstrahlung emission, the fractional variation in density should 
be half that of the surface brightness seen. This is not the case in 
reality, as projection effects are dominant. We constructed simple 
numerical simulations of 10 to 20 kpc wavelengt h density waves 
in a cluster density profile (we tried the profiles in lChurazov et al.l 
2003 and Sanders et al. 2004). The typical conversion factor from 
a fractional surface brightness to density perturbation is 2-3, with 
more suppression for smaller wavelength waves. 

We take the surface brightness image, and filter it with a high- 
pass filter as in Section [STI We use a looser filtering here as some 
of the observed 20 kpc waves are otherwise suppressed (not filter- 
ing anything below a wavelength 62 arcsec, increasingly filtering 
up to 124 arcsec and discarding everything longer wavelength). The 
original surface brightness ima ge was binned using the contour bin- 
ning algorithm ( ISandersll200 6) to have 10^ counts per region. We 
applied the same binning to the filtered image, and divided it by the 
binned surface brightness image. This created a fractional surface 
brightness variation map. The contour-binning routine follows the 
surface brightness closely, so bins are also aligned well with the 
ripples. 

For each pixel on the fractional surface brightness variation 
map, we compute the power in a spherical wave from Equation 
|2l assuming that the radius of the wave is the projected radius on 
the sky. Deprojected density and temperature values at each radius 
were c alculated from a fit to the average profiles in I Sanders et al.l 
1 I2OO4I) . In Fig. |4] we show a map of computed power values for 
each pixel, assuming a factor of 2 to convert from surface bright- 
ness variations to density variations. Many of the ripple features 
in Fig. [T] (right panel) are seen in this image, plus the radio bub- 
bles (which are not themselves sound waves). We bin the data in 
this image rather than use simple smoothing, as noise in the outer 
regions means that the power is overestimated. 
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Figure 6. Plot of mean power in sound wave and the cluster luminosity per unit length (in radius). The shaded region shows the the minimum and maximum 
sound wave power when the data are split into six equal sectors. The cluster luminosity shows the heating rate at each radius required to combat cooling. The 
vertical line shows the inner radius of our measurements due to the radio lobes. Note the difference in scales, the power being on the left in units of Iff*"* erg s^ ' 
and the luminosity on the right in units of lO''^ erg s^' kpc^' . 
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Figure 5. Plot of mean power in sound waves and the cluster cumulative X- 
ray luminosity inwards of where the cooling time is 4 or 5 Gyr. The vertical 
line shows the inner radius of our measurements due to the radio lobes. The 
total power in the lobes approaches lO'*' erg s"' jPunn & Fabianl2004h . 



In Fig.|5]we show the average values at each radius, masking 
out the central bubbles and north-west bubbles and the edges of the 
CCD to avoid filtering artifacts (using the regions in Fig.Q. On the 
plot we also display the cumulative luminosity from the cluster cal- 
culated fro m the deprojected den sity, temperature and abundance 
values from ISanders et alj i2004h . We accumulate the luminosity 
inwards from a radius of 75 or 100 kpc, the radii corresponding 



to wh ere the mean radiative cooling time of the gas jSanders et al.l 
is ~ 4 or 5 Gyr, the likely age of the cluster. Note that if the 
sound speed is a function of azimuth in the cluster (the temperature 
map indicates that this is so), then the phasing of the waves depends 
on azimuth. This will lead to smearing of the ripples in Fig.|5] 

We plot the wave power out to larger radii in Fig. |6] In this 
graph we also show the variation in wave power at each radius as 
a shaded region. This was computed by repeating the calculation 
of the average power in six equal sectors, and shading the region 
between the minimum and maximum values at each radius. At large 
radii only a couple of sectors were used, due to the position of the 
source on the detector. We show the luminosity of the cluster per 
unit length of radius in this plot to compare to the wave power. 

Fig. [5] shows that the net power in the ripples is a few times 
lO'*'* erg s^' and sufficient to offset a significant part of the radia- 
tive cooling within the innermost 70 kpc or so. The power implied 
by the analysis drops off with radius out to 120 kpc, with an e- 
folding length of about 50 kpc, consistent with models of viscous 
dissipation (Fabian et al. 2003a), which is required if heating by 
sound waves offsets radiative cooling. 

Larger power is seen near the edge of Fig.[6]around 105, 115 
and 125 kpc radius. These show possible evidence that the source 
was more powerful several 10^ yr ago. Such powerful shocks could 
have been created as several indiv idual bursts, or from a single burst 
producing multiple sound waves teriiggen et alj2007l) . We caution 
that most of the signal comes from a small angular region to the 
extreme SW of the main detector. Further observations covering a 
wider region are required to accurately determine the wave power 
at this radius. 

We also looked for structure in images of the cluster combin- 
ing all of the ACIS (Advanced CCD Imaging Spectrometer) CCDs. 
Fig. |7] shows an unsharp-masked image of the north of the clus- 
ter to large radii. We note that there appears to be a sharp edge at 
around 130 kpc radius in this direction (although it varies in radius 
in the northern sector). This is at approximately the same radius as 
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Figure 7. Further possible surface brightness discontinuities to the north of 
the cluster This is an unshaip-masked image, subtracting images smoothed 
with a Gaussian by 4 and 16 arcsec and dividing by the 16 arcsec map. 
Point sources were excluded from the smoothing. The apparent features are 
at radii of ~ 130 and 170 kpc (indicated by solid lines). 
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Figure 9. Metallicity map of the core of the cluster relative to solar Regions 
contain greater than 4 x 10** counts. Fits assume solar ratios of elements, 
but the results mainly depend on the iron abundance. Uncertainties for each 
spectral fit range smoothly from 0.025 in the centre to 0.07 Z,., to the ex- 
treme bottom left of this image. The circle shows the approximate position 
of the ancient bubble, and the diamond shows the high metallicity blob in 
its apparent wake. 



around 4 times larger than it was originally (if it retains its original 
energy). Rough scaling with the rising bubble NW of the nucleus 
suggests that it would have to have at least twice as much energy. 
The direction from the centre in which the dip is most visible (the 
longer solid line in Fi g. [7) is also directl y along the northern Ha 
filament and fountain jpabian et alj|200^ . If it is indeed a bubble 
then it shows that bubbles remain intact to very large radii in galaxy 
clusters. Bubble- like low pressure r egions were also seen to the 
south of the core jFabian et al .120061) . 

To constrain better the power in the feature at 130 kpc and 
confirm the radio bubble near 170 kpc radius, requires further deep 
observations by Chandra offset from the cluster core to improve 
the point spread function (the average radius of which is 10 arcsec 
at 10 arcmin off axis, compared with ~ 1 arcsec on axis). 



Figure 8. Fractional residuals from a /3 model fit to the surface brightness 
in sectors to the north at large radii. The full sector shows a profile between 
19 and 75° from west towards north, while the naiTow sector is between 59 
and 75°. 

the ripple seen in Fig. |6] The edge can be seen in residuals from 
/3 model plus background fits to a wide and narrow sector (Fig.[8l(. 
This discontinuity c an also be seen in an X MM-Newton observation 
of Perseus (Fig. 7 in lChurazov et ai]|2003h . 

At a radius of around 170 kpc is another feature. This appears 
to be a dip in surface brightness followed by a rise. This feature 
is particularly sharp in fits to the surface brightness in a narrow 
sector (Fig.[8). A possible interpretation is an ancient radio bubble 
which is still intact. The thermal gas pressure is likely to be around 
4 times less at this radius, so if the bubble remains intact it will be 



4 METALLICITY MAP 

Fig. [9] shows a metallicity map of the core of the cluster. This 
was gene rated by extracting spectra from contour-binned regions 
( ISandersll2006) containing ~ 4 x 10"* counts. The spectra were fit 
by a PHABS absorbed MEKAL model with the absorption, tempera- 
ture, emission-measure and metallicity free. Note that the plot does 
not clearly show the high abundance sh ell (possibly r narkin g the lo- 
cation of an ancient bubble) found by Sanders et alj feoOSi) as there 
were no new data in that region in these observations, and the spa- 
tial regions we use here are larger. 

To examine the variation more quantitatively, we have re- 
peated the spectral fitting using regions containing ~ 2.5 x 10^ 
counts to decrease the size of the uncertainties. The radial temper- 
ature and metallicity variation are plotted in Fig.[TOl generated by 
plotting the average radius of each bin against the value obtained 
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Figure 10. Radius versus projected temperature and metallicity. Values 
were measured from spectral fitting a MEKAL model to spectra from bins 
with ~ 2.5 X 10^ counts. Also plotted is the projected entropy as a function 
of radius (see Section l4~l1 . 

from that region. At each radius there is a large spread in tempera- 
ture and abundance. 

4.1 Metallicity relations 

In ISanders et al.l ( |2004) we plotted the temperature of regions 
against their metallicity. As the temperature of the gas declines to- 
wards the centre, the metallicity reaches a maximum at a radius of 
40 kpc, then decreases again. We plot the metallicity-temperature 
relation from the new data in Fig.[TT](top panel). Although we see 
a similar relation, the reduced errors bars show that there is sig- 
nificant metallicity scatter at each temperature. Indeed the scatter 
appears similar to that in the radial plot of the metallicity (Fig.llOt. 
The temperature of the gas does not appear to correlate better with 
the metallicity than the radius. 

Another interesting physical quantity is the entropy of the gas. 

-2/3 

Entropy in clusters is usually defined as A" = kTiig . Using the 
XSPEC normalisatior[3 per unit area on the sky. A' <x n^d, where d 

' XSPEC MEKAL and APEC normalisations are defined as 
{10"''*/«enHdV}/{47r[DA(l + z)]^}, where the source is at redshift 





Projected entropy 



Figure 11. Temperature versus metallicity and pseudo-entropy versus 
metallicity plots. All values are projected. The bins used are the same as 
in Fig.fTol 

is a depth in the cluster, which we assume to be constant, we can 
calculate a pseduo-projected entropy quantity A' The plot 

of this quantity against the metallicity (Fig.[TT]bottom panel) looks 
similar to the temperature-metallicity plot, with a great deal of scat- 
ter at each entropy value. The values use A' in XSPEC normalisation 
units per square arcsec and kT in keV. 

The scatter in the abundance seems unrelated to the tempera- 
ture, radius, or entropy. Higher metallicity regions will have shorter 
mean radiative cooling times as the line emission is stronger, but 
this effect is not strong at temperatures of ~ 3 — 7 keV, so it is per- 
haps not surprising that metallicity and temperature or entropy (and 
therefore cooling time) are unrelated. 



4.2 Central abundance drop 

There is also a drop in metallicity in the central regions (as can be 
seen in the metallicity map in Fig.|9). This feature appears to be un- 
affected by including other components, such as extra temperature 
components, or powerlaw models. 

We have investigated whether the central drop in metallicity 
could be caused by projection effects. We use the spectral depro- 
jection method outlined in Appendix lAlto construct a set of depro- 
jected spectra. Fig. [12] shows the iron metallicity profile computed 



z and angular diameter distance D^, and the electron number density nc 
and Hydrogen number density nn are integrated over volume V. 
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Figure 12. Iron metallicity profile to the NW of the cluster between position 
angles 290.0 and 4.7° from 3C 84. This was produced by fitting VMEKAL 
models to projected and deprojected spectra. 
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by fitting isothermal VMEKAL models to the spectra before and af- 
ter deprojection. In the fit the elemental abundances of O, Ne, Mg, 
Si, S, Ar, Ca, Fe and Ni were allowed to vary, with the gas tem- 
perature, the model emission measure, and the absorbing column 
density. The analysis shows the peak at around 40 kpc radius is en- 
hanced by accounting for projection and the central drop remains. 

The drop in metallicity in the central regions appears to be 
robust. The entropy of the gas in the central regions is much lower 
than further out, which means it is difficult for the source of the low 
metallicity gas to be from larger radius. 

4.3 Small scale variation 

Much of the structure in the metallicity map is real. To demon- 
strate this, we show in Fig.[T3]a metallicity and temperature profile 
across the high metallicity blob apparently in the wake of a pos- 
sible ancient bubble down to the extreme SW (marked on Fig.|9] 
ISanders et al. ||2005). The metallicity profile shows a significant 
peak with a width of '--^ 5 kpc at the position of the blob (we also see 
~ 5 kpc metallicity features near the filaments in Section[6]l. There 
is no obvious correlation between the temperature and metallicity 
profiles. Taking the global diffusio n coefficient of 2 x 10^^ cm^ s^' 
measured bv Rebusco et al.l ( l2005h in Perseus and a lengthscale of 5 
kpc, the lifetime of such a feature is only ~ 40 Myr. This is roughly 
comparable to the buoyancy timescale estimated for the ancient 
bubble of 100 Myr (Dunn et al. 2006), given the large uncertain- 
ties, but may require that diffusion is suppressed on small scales as 
it appears unlikely that the metals could have been injected in-situ. 
The sharp edges of the feature would imply diffusion times of only 
a few Myr, requiring significant suppression. 

The blob has a metallicity approximately Q.lXr., higher than 
its surroundings (which are at approximately 0.6 Zp). This is prob- 
ably a lower limit because of projection effects. Assuming the blob 
is sphere with diameter 5 kpc, and a local electron density of around 
0.037 cm^-' ( Sande rs et al.ll2004h . it represents an enhancement of 
around 2.6 x 10^ Mq of Fe. If most of this enrichment is due to 
Type la supernova, then this corresponds to around 3.7 x 10'* su- 
pernova (assuming 0.7 Mp of Fe produced per supernova). Taking 
the timescale from diffusion, this corresponds to 0.1 supernova type 



Figure 13. Abundance and temperature profiles across the high metallicity 
blob near 3!"19!"40?, +41°29'15". The profile crosses the blob in the SE to 
NW direction, inclined by an angle of 40° from the west northwards. 



la per century for this 5 kpc radius region, which is a significant 
fraction of the rate expected from a single galaxy. 



5 THE HIGH VELOCITY SYSTEM 

iGillmon et al] ( |2004|) studied the High Velocity System (HVS), a 
distinct emission-line system at a higher velocity than NGC 1275 
(iMirtkowski 195Z), using an earlier 200-ks Chandra observation 
of the system. Th e system is observed in absorption in X-rays 
dFabian et al.l200Cl), as it lies between most or all of the cluster and 
the observer. iGillmon et ai] ( |2004|) mapped the absorbing column 
density and placed a lower limit of 57 kpc of the distance between 
the nucleus of NGC 1275 and the HVS. They obtained a total ab- 
sorbing gas mass of 1.3 X 10''M,;„ assuming Solar metallicities. 

Here we repeat the mapping analysis with the new 900-ks 
dataset. Fig.[T4](top panel) shows an image in the 0.3-0.7 keV band, 
clearly showing the absorbing material. To quantify the amount of 
material we fitted an APEC thermal spectral model for the cluster 
emission, absorbed by a PHABS photoelectric absorber to spectra 
extracted from a grid of 0.74 x 0.74 arcsec regions around the HVS. 
In the fits the temperature of the thermal emission, its normalisa- 
tion, and the absorbing column density were free. The abundance 
of the thermal gas was fixed to 0.7 Zq (based on the results of spec- 
tral fitting to larger regions near the core of the cluster; the results 
are very similar if this is allowed to be free). Spectra were binned 
to have a minimum of 5 counts per spectral channel. We minimised 
the C-statistic to find the best fitting parameters (note that no back- 
ground or correction spectra were used in the fits here, as the cluster 
is much brighter than the background here). 

In Fig. [14] (bottom panel) we show the resulting Hydrogen 
column density map over the HVS region, including the Galac- 
tic contribution. We stress that these values are equivalent Hy- 
drogen column density. The measurements are most sensitive to 
the abundance of Oxygen in the absorber. The Hydrogen col- 
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Figure 14. (Top panel) Image of the high- velocity system in the 0.3 to 
0.7 keV band smoothed by a Gaussian of 0.49 arcsec. (Bottom panel) 
Measurements of the equivalent Hydrogen column density over the high- 
velocity system region, in units of 10^~ cm"~ (subtracting Galactic absorp- 
tion of 0.12 X 10"^ cm^^). The values were measured from spectral fitting 
in 0.74 X 0.74 arcsec pixels. The central nucleus of 3C 84 shows up as a 
point-like object near 3!'19"'48f, +4r30'40". The uncertainty on the col- 
umn density on individual pixels is around 0. 1 for the regions of highest 
absorption and 0.03 outside of the HVS. 

umn density is calculated assuming the Solar abundance ratios of 
lAnders & Grevessd jl989h . which gives an Oxygen to Hydrogen 
number density ratio of 8.51 x 10^^. The column density is flat 
in this region beyond the HVS. We obtain an average value of this 
Galactic component of 0.12 x 10^^ cm^^ , consistently from several 
surrounding regions. 

The total number of absorbing atoms can be calculated us- 
ing the mean absorption over the HVS, the Galactic contribution, 
and the distance to Perseus. We converted this into a total ab- 
sorbing mass of 1.1 X 10'' Me, assuming Solar metallicity ratios. 
This is lower than the value of (1.32 ±0.05) x 10^ Mg quoted by 
iGillmon et al.l . We repeated our analysis using the 1.96 arcsec bin- 
ning factor used in that paper, but this did not significantly change 
our result. The twenty per cent difference may be due to the differ- 
ent calibration used in the earlier analysis, particularly the lack of 
correction for the variation in the contaminant on the detector, as 
the Galactic value is important to determine the total mass. 

To place an improved lower limit of the distance of the HVS 
from the cluster nucleus, we examined a 0.3 to 0.7 keV image, 
where the absorption is strongest. Taking 1-1.5 arcsec diameter re- 
gions (with total area 3.9 arcsec^) in the three highest absorption re- 
gions we find 17 counts in total in this band. Exposure-map correct- 
ing this surface brightness, and comparing it to the exposure-map 
corrected image of the cluster, we find that the surface brightness 
from the cluster only goes down to this level at a radius of approx- 
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Figure 15. Location of the regions used to measure the thermal propeilies 
across the filaments, displayed on a 0.3 to 1 keV image. Each of the 80 
regions is 1 x 27.3 arcsec in size. 

imately 110 kpc. This is a lower limit of the distance of the HVS 
from the cluster centre, as if the HVS were closer, the cluster emis- 
sion should be 'filling in' the de crement in count rate caused by the 
absorption jGillmon et alJl20()4l) . This lower limit is a significantly 
better the previous value of 57 kpc. 

As we discuss in Section 17.1.1 1 the HVS could have a shock 
cone behind it if it is travelling exactly along our line of sight to- 
wards the cluster. If this is the case, the shocked material would 
have a significantly lower density than the cluster at the same ra- 
dius. This would mean that the lower limit we compute above 
would be overestimated. 



6 PROFILES ACROSS FILAMENTS 

To investigate the thermal structure of the X-ray filaments associ- 
ated with the Ha nebula in detail, we extracted spectra from small 
1 X 27.3 arcsec boxes in a profile across the filament (the regions 
we used are shown in Fig.|15t. We created background spectra, re- 
sponses, ancillary responses, and out of time background spectra 
for each region. 

6.1 Multiphase model 

Our first model was to fit the spectra with a multiphase model con- 
sisting of APEC thermal components at fixed temperatures of 0.5, 
1, 2, 4 and 8 keV. The normalisations were allowed to vary and the 
metallicities of each of the components were tied together. They 
were absorbed with a PHABS photoelectric absorber which was al- 
lowed to vary. We assume each component has the same metallic- 
ity as we cannot measure them independently. The measurement 
is likely to be driven by the cooler components as they are line- 
domi nated. The model i s similar to that used to produce figure 
12 in I pabian et al.l | |2006|) , mapping the multiphase gas. We show 
the emission measure of each temperature component in Fig. \16\ 
with the absorbing column density and metallicity. Also shown is 
an unsharp-masked image of the region examined in the 0.5-7 keV 
X-ray band (the X-ray image has been rotated so that the 1 -arcsec 
bins lie horizontally across the region) and the Ha image from 
IConselice etal.l ( l200ll) taken using the WIYN telescope. We chose 
the rotation angle so that the bins were lined up across the central 
X-ray features. The Ha filaments are not quite aligned in angle. 
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The plot shows that the filaments do not contain gas at temper- 
atures of ~ 4 keV. At 2 keV we start to see the filaments, although 
their signal is fairly weak against that from nearby gas. The fila- 
ments are very strong near 1 keV, and still visible at 0.5 keV. There 
are some point-to-point differences however. The strong filament 
at 39 arcsec distance is strong in 1 and 0.5 keV, but the one nearby 
46 arcsec does not appear at 1 keV, but does at 2 and 0.5 keV. 

The column density profile is fairly flat. A linear model fitted 
to the column density profile gives a reduced of 86/78 — 1.10. 
We see no evidence for additional X-ray absorption associated with 
the filaments. 

The metallicity profile shows several features, which matches 
the complex structure seen in the metallicity map (Fig. |9}. A sim- 
ple linear fit gives a reduced of 105/78 = 1.34. There are quite 
strong peaks at the 59 arcsec, and at 12 arcsec. Neither is exactly 
at the position of a filament, though the one at 59 arcsec is offset 
by an arcsec or two from some gas from 2 to 0.5 keV. It is possible 
that some of the metallicity variation is caused by incomplete mod- 
elling of the multiphase gas within the filaments, but we observe 
almost identical variation with single phase and cooling flow (be- 
low) models. All models assume that the gas at each temperature 
has the same metallicity, though the cooler gas will dominate in the 
measurement of the metallicity. 

6.2 Cooling flow model 

Given the range in the temperature in the filaments, we try fitting a 
cooling flow model to the spectra, as the gas may be cooling within 
the filament. We used the isobaric MKCFLOW model, which models 
cooling between two temperatures at a certain metallicity and gives 
the normalisation as a mass deposition rate, plus a MEKAL single 
phase thermal component. The upper temperature of the cooling 
flow component was tied to the temperature of the thermal compo- 
nent, and the metallicities were also tied. Both components were 
absorbed by a PHABS absorber. Fig.[T7]shows the mass deposition 
rates obtained with fixed minimum temperatures. The 0.25 keV rate 
is equivalent to a traditional cooling flow where the gas cools out of 
the X-ray band, while the others are truncated cooling flows. Also 
plotted is the best fitting minimum temperature if it is allowed to 
vary. 

The results show that within the filaments, the model adds 
increasing amounts of the cooling flow component, indicating a 
range of temperatures. Using a cooling flow component with a 
0.25 keV minimum temperature, the emission measure of the sin- 
gle phase thermal component varies relatively smoothly over the 
region. There is less apparent co oling as the minimum temperature 
decreases (as with the results of IPeterson et alj|2003h . until below 
~ 0.8 keV where the results are approximately consistent. If the 
gas is actually cooling in these filaments, these results indicate a 
rate of '--^ 0.25 M,;, yr^' in the central filament in this small region 
examined. All of the ~ 0.5 keV X-ray emitting gas observed in the 
Perseus cluster appears to be as sociated with the H o: emitting fila- 
ments or is close to the nucleus jpabian et alj2006t) . The lower gas 
temperature measured here is consistent with zero in the centres 
of the filaments, but appears to increase by a few arcsec. We note 
that the best fitting absorbing column density is almost identical to 
the values from the multiphase model in Fig. [T6l so no additional 
absorption is required. 

In Fig. [TS] we compare the flux from the cooling flow 
component cooling down to 0.25 keV (above) against the flux 
in the continuum-sub tracted Ha waveband from the data from 
IConselice et al.l ( 1200 1[) (using the fluxes given in several regions to 



calibrate the count to flux scale). The Ha filter used had a central 
wavelength of 6690 A and a FWHM of 77 A, and the fluxes have 
been corrected for Galactic extinction, and includes N II emission. 
The plot shows that the X-ray flux of the filaments, assuming a 
cooling flow model, is very close to the Ha-l-N [II] flux for most of 
the filaments. The flux from the multitemperature model (Fig.|16ll 
is somewhat similar to that from the cooling flow model. If the 2, 4 
and 8 keV components are ignored, the flux is smaller by a factor of 
~ 4 than the cooling flow flux, but including the 2 keV flux boosts 
this to well above this value. 



6.3 Filament geometry 

If the gas in a filament is in pressure equilibrium with its surround- 
ings, the pressure is known and thus is the temperature of the gas in 
the filament, so we can estimate the volume of the emitting region. 

We fitted a simple absorbed two-temperature MEKAL spectral 
model to the 1 x 27.3 arcsec bin at 38 arcsec offset where the fil- 
aments are strongest (see Fig. [T6t . The best fitting temperatures 
from the spectral fit are 0.69 ±0.05 and 3.36 ±0.09 keV, corre- 
sponding to the temperature of the filament and its surroundings, 
respectively, under the assumption that the gas has a single temper- 
ature in the filament. Taking the emission measure of the filament 
component and an electron pressure of 0.128keVcm~-^ (from the 
deprojected values in Sanders et al i 20041 . Fig. 19), then the volume 
of the emitting region is 1.0 X lO*'^ cm . The extracted region on 
the sky is 0.37 x 10.2 kpc. If we assume the filament is a cuboid 
with the facing side of the dimensions given, we calculate a depth 
of 9 pc. If the filament is a single long and thin tube with length 
10.2 kpc, it would have dimensions of ~ 60 pc across. 

The lack of depth shows the filament is unlikely to be a sheet 
viewed from the side, and the disparity with the other dimensions 
suggests that the filament is in fact made up of small unresolved 
knots of X-ray emission. If there were lO"* such blobs, they would 
have dimensions of ~ 30 pc. 



6.4 Flux emitted from the filament 

Taking the volume above for the filament in the 38 arcsec bin and 
assuming pressure equilibrium implies there are ~ 1.9 X lO''^ elec- 
trons in the cooler filamentary component. The difference in tem- 
perature of the cooling X-ray emitting gas from the surrounding gas 
is ^ 2.7 keV. This implies there would be ~ 2.5 x lO^'* erg released 
if the X-ray emitting filament was to cool out of the surrounding in- 
tracluster medium (assuming twice the number of particles as elec- 
trons and 3/2 fcr energy per particle, which does not include any 
work done by the surroundings on the cooling filament). 

The luminosity of the bin in the Ha waveband is 6 x 
10"*" erg s"' .This can be multiplied by a factor of 20 to account for 
other line emission, for example Lya. Therefore if the line emission 
is powered by cooling out of the intracluster medium, the timescale 
for this process is ~ 7 x 10** yr. The dynamical timescale for a 30 pc 
region, assuming a sound speed of 300 km s^' , would be ~ 10^ yr, 
so it is approximately possible for the cooling medium to be replen- 
ished. 

Conduction of heat from the surrounding intracluster medium 
may also be able to power the filaments. For a 2 keV plasma, the 
Spitzer ( 1962) conductivity is ~ 9 x lO" ergs^' cm" we 
assume a geometry of a 30 pc wide tube which is 10.2 kpc long and 
if the heat is being conducted between 3.4 and 0.7 keV over 30 pc, 
this corresponds to a heat flux of '-^ 5 x 10^^ erg s^' . Such a flux is 
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sufficient to fuel the line emission. It however depends greatly on 
the assumed conductiv ity (which depends on temperature to the 5/2 
powe r) and geometry jBohringer & Fabiarilll989l : lNipoti & BinnevI 
|2004 . To provide the energy for the line emission, the heat must be 
able to travel to presumably smaller and cooler regions than that of 
the the 0.7 keV gas, which conductivity makes it increasingly hard 
to do. Note that any conduction model requires that the surround- 
ing soft X-ray emitting gas and the filament have a low relative 
velocity. If conduction is too efficient at transporting heat the fila- 
ment wi ll evaporate, leading to a c ertain critical minimum size for 
growth teohringer & Fabianlll989h . which depends on how much 
conduction is suppressed and the temperature and density of the 
surrounding ICM. 



7 HARD X-RAY EMISSION 

As discussed in the introduction, ISanders et all ( l2004b found evi- 
dence for a distributed hard component surrounding the core of the 
cluster by fitting a multitemperature model with a 16 keV com- 
ponent. The total 2-10 k eV luminosity of thi s hard component is 
~ 5 X 10"*^ erg s^' . Later lSanders et al. 1(2005") used a thermal plus 
powerlaw model to fit the data giving a similar flux. This hard com- 
ponent has been confirmed with XMM-Newton data (Molendi, pri- 
vate communication). Assuming this emission is the result of in- 
verse Compton scattering of IR and CMB photons by the popula- 
tion of electrons which also emit the observed radio, the magnetic 
field over the core of the cluster was mapped. 

Here we examine the hard flux from the cluster using this 
very deep set of observations, and investigate the effect of differ- 
ent spectral models. We binned the data using the contour binning 
algorithm with a signal to noise of 500 (~ 2.5 x 10^ counts in each 
bin). 

Our first model is a multitemperature model made of differ- 
ent temperature components plus a powerlaw to account for the 
hard emission. This is a more complex model than that used by 
ISanders et all ( |2005|) , as we wished to account for the known cool 
gas in the cluster which may affect the powerlaw signal if not mod- 
elled correctly. Hotter gas projected from larger radii in the cluster 
can also give a false signal. The data were fitted using a model 
made up of APEC thermal components at fixed temperatures of 0.5, 
I, 2, 3, 4 and 8 keV, plus a F = 2 powerlaw, all absorbed by a 
PHABS photoelectric absorber. The normalisations of the compo- 
nents were allowed to vary, the metallicities tied to the same free 
value, and the absorption was free. We used a F = 2 powerlaw here 
as this is close to the best fitting photon index of the radio emission 
I Siibrin jl993h, and the best fitting X-ray spectral index in the core 
I Sanders et al. l2005h . Fig. [19] shows the emission measure per unit 
area for each of the thermal components, the powerlaw normalisa- 
tion per unit area, and the broad band Chandra image of the same 
region. 

The thermal gas maps show similar distributions to the 
earlier analyses with smaller uncertainties (" Sanders et J] |2004 
lFabianetalj|2006l) . We plot the radial profile of the 2-10 keV 
flux of the powerlaw component in Fig. [20] Also shown is the ra- 
dial profile for a F = 1.5 powerlaw instead of the F = 2 power- 
law, or a 16 keV th ermal component, and the older results from 
ISanders et al.l ( |2005| ) (after subtracting the estimated 'background' 
from hot projected thermal gas). The total emitted flux from each 
of the models is fairly consistent, except at large radii where the 
signal is low. Hot thermal gas, as expected, gives similar results to 
a F = 1.5 powerlaw. The best fitting powerlaw index looks similar 



to the results in lSanders et al.l j2005l) , with a transition of F ~ 2.2 
powerlaw in the centre to F ~ 1 .4 in the outskirts. 

We map the distribution of hard flux per unit area for a vari- 
ety of different models in Fig. [2T] We show the variation of flux 
(top left) just using a single thermal model plus powerlaw (as in 
ISanders et al.ll200a) , (top right) fitting that model just in the high 
energy band, and (second row left) allowing the powerlaw index to 
vary. We also show the results (second row right, and third row) 
from multitemperature plus powerlaw models with F = 2, 1.8 and 
1.5, and (bottom left) fitting for the index. Finally we show the re- 
sult (bottom right) using a multitemperature model including a 16 
keV hard component. 

There are problems with steep powerlaw components as they 
predict significant extra flux at low X-ray energies. Such flux is not 
observed, and therefore the best fitting photoelectric absorption in- 
creases in the central regions where the powerlaw is strong. It is 
possible that there is excess absorption in the central regions, but 
the flux of the powerlaw is very closely correlated with the absorp- 
tion, so some physical connection between the two would be re- 
quired. Presumably the non-thermal emission process would need 
to be dependent on absorbing material, which appears unlikely. 

Fig. [22] shows the absorption distribution for different mod- 
els. The top-left panel shows the distribution from fitting a single 
MEKAL plasma to the projected spectra. There is obvious variation 
across the image. Some of this variation may be due to the buildup 
of contaminant on the ACIS detector, but the current calibration 
should account for this in the creation of ancillary response matri- 
ces. Probably most of the variation is because the cluster lies close 
to the Galactic plane (i '-^ — 13°). If the image is aligned to Galactic 
coordinates, the variation is mostly in Galactic latitude. 

Fitting using a multitemperature model plus a F = 1.5 power- 
law or 16 keV thermal component produces absorption maps very 
similar to the single temperature map. These models require no ad- 
ditional absorption. The F = 1.8 requires moderate additional ab- 
sorption, and F = 2 produces absorption clearly correlated to the 
powerlaw flux (Fig.[2T]centre-top row, right column). 

To demonstrate the effect of the hard component on the spec- 
tral fit, we show in Fig.[23]the contribution of the hard component 
to the best fitting model to a spectrum from a region around 1.8 ar- 
cmin north of the nucleus. For a 16 keV component the effect is 
around 10 per cent at low energies, increasing to fifty per cent at 
high energies: the hard component is about one third or more flux 
above 4 keV. 



7.1 Origin of the hard emission 

If we assume that the hard component is real and not an instrumen- 
tal or modelling artifact, there are two sets of possible emission 
mechanisms, thermal or non-thermal. If it is non-thermal emis- 
sion it is likely to be inverse Compton emission from CMB or 
IR photons being scattered by relativistic electrons in the ICM 
(Sanders et al. 2005). Another possible origin is hot thermal gas. 
The probable origin of such material in the cluster is from a shock, 
and the obvious candidate is the HVS. 



7.1.1 Thermal origin: merger of HVS 

From Section[5] the HVS lies at minimum distance from the clus- 
ter core of 110 k pc, where the electr on density of the gas is 
~ 5 X lO^^cm"^ dSanders et J]|2004l) . The HVS is moving at 
3000 km s^' relative to the main NGC 1275 system along the line 
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Figure 19. Results from the multitemperature plus F = 2 powerlaw fits to regions containing ~ 2.5 y 
unit area for each of the components, and finally the broad band Chandra image of the same region. 
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of sight to the observer jMinkowskil 1 1957|) . Assuming a 5 keV 
plasma for the cluster, and that the galaxy is travelling along the 
line of sight, the Mach number of the HVS is around 2.3. 

The collision should produce a shocked Mach cone around the 
HVS. If the HVS lies at its minimum distance from the cluster, us- 
ing the Rankine-Hugoniot jump conditions, the density of the post- 
shocked material should be a factor of 2.1 times greater than its 
surroundings {rie ~ 10^^ cm^^ ) and its temperature around 15 keV. 

If the origin of the hard component is the shocked material 
in the cone, and it lies at the minimum distance from the cluster, 
we can estimate the depth of the cone from the emission measure 
of the 16 keV component (close to the expected 15 keV tempera- 
ture) using the density above. The depth we derive from the peak 
of the emission is between 200-300 kpc. If this is the correct in- 
terpretation, the layer of shocked material is extremely large. This 
number can be reduced significantly if the preshocked material has 
a higher density, as the depth is inversely proportional to the square 
of the density. At the previous minimum distance of around 60 kpc 
from the HVS from the core of the cluster dGillmon et alj|2004h . 
the minimum depth is only a few kpc, as the density was increased 
by a factor of 10. If the HVS lies further away than our minimum 
distance from the cluster core, it becomes much harder for the ther- 
mal interpretation to be a plausible explanation. Whether a thermal 
origin for the hard component is plausible depends on whether our 
lower limit of the distance from the cluster to the HVS is overesti- 
mated (see Section|5]l. 



7.1.2 Non-thermal origin: inverse Compton processes 

If the origin of the flux is from inverse Compton emission, there 
is a large population of relativistic electrons scattering CMB or IR 
photons. The bright radio emission from the mini radiohalo in the 
cluster core indicates that there are relativistic electrons, but we 
cannot directly observe the required y ~ 1000 electrons by their 
synchrotron emission. The powerlaw index of the inverse Compton 



emission will be the same as the radio emission, if a single popu- 
lation of electrons produces both. We do not know whether there is 
a single population. There could be a break in the powerlaw index, 
or a cut-off in the population. 

Under the assumption of a single population, we previously 
estim ated the magnetic fi eld in the core of the cluster to be a 0.3- 
IjuG l lSanders et all2005h . 

There are some potential issues with an inverse Compton ex- 
planation. If the source electron population creates a F ~ 2 power- 
law, we require significant amounts of additional absorption in the 
core of the cluster (Fig. l22t . as the powerlaw is significant at low 
X-ray energies, unless it breaks in the X-ray band. Such absorp- 
tion is not seen on the X-ray spectrum of the nucleus using XMM 
dChurazov et alj2003h . The excess absorption is not required with a 
flatter powerlaw index (F ~ 1.5). However a flat F = 1 .5 powerlaw 
will emit significant flux in the hard (> 20 keV) band. 

Hard X-rays w ere observed from Perseus using HEAO 1 
jPrimini et al.lll98l[) in the 20-50 keV band with a photon index of 
around 1.9. The flux translates to a total luminosity in the 2-10 keV 
band of around 1 .6 x 10^^ erg s^' above the thermal emission. This 
component was f ound to be variable over a four year timescale. 
iNevalainen et al.l ( l2004h observed a flux around four times lower 
than this value using the high energy PDS detector on BeppoSAX. 
The variable component must be associated with the central nu- 
cleus rather than inverse Compton emission. The separation of the 
nuclear spectrum from any hard compone nt is difficult, particularly 
if the nucleus varies on short timescales. INevalainen et al.l ( l2004h 
concluded that the central nucleus can account for all of the non- 
thermal emission observed from Perseus using BeppoSAX. More 
sensitive measurements of the hard flux from, for example, Suzaku 
are vital to resolve this issue. Observations with higher spatial res- 
olution close in time would help remove uncertainty about the nu- 
clear component. 

The cooling time for electrons producing 10 keV electrons 
from inverse Compton scattering of CMB photons is '--^ 10*^ yr It 
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is therefore likely that the spectrum could be broken at these en- 
ergies. Such breaks could reduce the need for increased absorption 
for steeper powerlaw models here, or reduce the hard X-ray flux for 
those with flatter spectra. 

If the emission is the result of inverse Compton emission, the 
pressure of the relativistic electrons, P, is rel ated to the emissiv ity 
of the inverse Compton emission, S, by (e.g. lErlundetai]|2006l) 



where U is the energy density of the photon field being scattered, 
y is the Lorentz factor of the electron scattering the photon to the 
observed waveband, nig is the rest mass of the electron. Ox is the 
Thomson cross-section and c is the speed of light in a vacuum. If 
we assume that the depth of an emitting region in the hard X-ray 
maps is its radius, and that th e X-ray emission is the result of scat- 
tering CMB and IR photons dSanders et aDl2005h we can estimate 
the nonthermal pressure as a function of radius. We plot this in 
Fig- US showing that the nonthermal pressure is comparable to the 
thermal pressure near the centre of the cluster. 



8 THERMAL CONTENT OF THE RADIO BUBBLES 

Earlier Chandra data have been used to limit the amount of therma l 
material material within the radio bubbles ( ISchmidtetal.ll2002h . 
This analysis is made more difficult because the geometry of the 
core of the cluster is complex, and techniques accounting for pro- 

i' ection appear not to work generally when examining the bubbles 
Sanders et alj2004h . Here we place stringent limits on the volume 
filling factor of thermal gas using this 900-ks combined dataset us- 
ing a comparative technique which depends far less on geometry. 

We fit the projected spectrum from the inside of the bubble 
with a model made up of multiple temperature components at rel- 
atively low temperatures (fixed to 0.5, 1, 2, 3, 4 and 6 keV with 
normalisations varied and metallicities tied together) to account for 
the projected gas, plus a component fixed to a higher temperature 
to test for the existence of hot thermal gas within the bubble. We 
also do the fit again with an additional F = 2 powerlaw to account 
for any possible non-thermal emission. We compare the normali- 
sation per unit area of the component accounting for hot thermal 
gas in the bubble against that from a fit to a neighbouring region 
at the same radius in the cluster. The temperature of the gas in the 
bubble is stepped over a range of temperatures in the bubble and 
comparison regions. 

The normalisation per unit area in the bubble and comparison 
region can be converted into an upper limit on the difference in nor- 
malisation per unit area between the two. We use any positive dif- 
ference between the bubble region normalisation and background, 
plus twice the uncertainty on the difference (Using the positive un- 
certainty on the background and the negative uncertainty on the 
foreground. This is slightly pessimistic compared to symmetrising 
the errors), to make a 2a upper limit. 

Assuming a volume for the bubble, an upper limit on the den- 
sity of gas at that temperature can be calculated, assuming the gas 
is volume filling. If the gas is at pressure equilibrium with its sur- 
roundings and the pressure is known, then a limit of the volume 
filling fraction of gas at that temperature can instead be calculated. 

We examine the inner SW bubble and ghost NW bubbles using 
the regions shown in Fig. |25] Also indicated are the regions used 
for background. The inner NE bubble is obscured by the High Ve- 
locity System, so we do not consider that. The ghost S bubble has 



somewhat uncertain geometry. We try to place the regions away 
from any contamination by low temperature gas (though we have 
tried alternative regions with little effect), and the background re- 
gions at a similar radius to the bubbles. We assume the bubble re- 
gions are cylindrical in shape, with depths of 9.4 kpc and 13.4 kpc 
for the inner SW and ghost NW bubbles, respectively. We take 
the deprojected electron pressure of the surrounding thermal gas 
from the mean value of the sec tors surrounding the bubbles in fig- 
ure 19 of Sanders et al. Il l2004h . This leads to values of 0.195 and 
0.111 keV cm^^ for the inner SW and ghost NW bubbles, respec- 
tively. 

For the two bubbles, including and not including the F = 2 
powerlaw component, we show the normalisation (emission mea- 
sure) per unit area for the foreground and background regions in 
Fig. [26] We convert these values to upper limits on the volume- 
filling fraction in Fig. [27] We find that if there is volume-filling 
thermal gas within the bubbles, it must have a temperature of less 
than '--^100 keV for either bubble. 



9 DISCUSSION 

We have investigated several issues from the Chandra data of the 
Perseus cluster and now attempt to tie the interpretation of the vari- 
ous phenomena more closely together. In particular we try to under- 
stand the interplay between heating and cooling in the cluster core, 
and how energy is transported and distributed through the ICM. 

9.1 Sound waves 

Observations of Perseus and many other similar clusters show that 
jets feed relativistic plasma from around the massive central black 
hole into twin radio lobes. The power associated with this process 
is large and comparable to that required to balance radiative cool- 
ing within the regions of 50-100 kpc where the radiative cooling 
times is 3-5 Gyr (Rafferty et al. 2006; Dunn & Fabian 2006). In the 
case of the Perseus cluster the jets are producing on average around 
5 X lO'*^ — 10^^ ergs^' in total, as estimated from the PAV work 
done over the 5-10 M yr age of the bubbles ( iFabian et alj|2002bl ; 
jPunn & Fab"iarill2004h . These estimates depend on the filling fac- 
tor of the bubbles by the relativistic plasma, which from the results 
given in Section[8]could well be unity. Therefore at least about 20- 
40 per cent of the bubble power goes into sound waves. 

We expect the jets to provide power more or less continuously 
over the lifetime of the cluster core (several Gyr). This is supported 
by the hi gh incidence of bubble s found in cluster cores that require 
heating ( jPurm & Fabian] I2OO6I) and by the train o f ghost bubbles 
seen in our Perseus cluster data jFabian et al.l200d l). Note however, 
that there is probably much fluctuation in the jet power over short 
timescales (the radio source has weakened in strength over the past 
40 yr), but variations over a cooling time of 10** yr are probably less 
than an order of magnitude. 

We now address the question of how the energy in the bub- 
bling process, occurring on a timescale of 10 Myr, is fed into the 
bulk of the core and dissipated as heat which balances the cool- 
ing. The smooth cooling time profiles seen in cluster cores, and the 
peaked metal distributions argue for a relatively gentle, continuous 
(on timescales of 10^ — 10^ yr or more) and distributed heat source. 
The inflation of the bubbles does PdV work on the surroundings 
and thus creates pressure waves - sound waves - which carry the 
energy outward in a roughly isotropic manner. We discovered rip- 
ples in surface brightness in the first 200 ks of the Perseus cluster 
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data jpabian et alj|20033) which we interpreted sound waves and 
have extended the analysis of them here. 

Fig. m shows that there is considerable power in these waves, 
around 2 — 3 x 10^"* erg s^' at radii of 30-70 kpc. It is similar to 
the level of heat required to offset radiative cooling within that re- 
gion. The wave power drops with radius between 30 and 100 kpc 
indicating that the energy is being dissipated and confirming the 
idea that viscous dissipation of the sound waves is the distributed 
heat source. The dissipation length is comparable to that estimated 
on the basis of Spitzer-Braginski viscosity (Fabian et a l. 2003a), al- 
though the process in the likely tangled magnetic, and cosmic-ray 
infused, plasma may be more complex and involve some form of 
bulk viscosity. 

There is however a deficit of power shown by this analysis 
within 30 kpc. Although it could just be variability of the central 
power source (the jets), our results on the nonthermal component 
associated with the radio minihalo indicates that there is a signif- 
icant nonthermal pressure from the cosmic rays there, comparable 
to the thermal pressure. This would quadruple (if the nonthermal 
wave pressure is the same as the thermal pressure in Equation |2]l 
the predicted wave power in this region, bringing it into agreement 
with the region between 30 and 60 kpc and meaning that dissipation 
of sound waves can be the dominant heat source balancing radiative 
cooling. 

The increase in power seen in Fig.[6]beyond 100 kpc can be ex- 
plained as associated with the power ~ 10** yr ago being 2-3 times 
larger than present. This may fit in with the pre sence of the long 
Northern optical filament jConselice et al. 2001 ), which if drawn 
out from the centre by a bubble i Hatch et alj2006l) must have been 
an exceptional bubble, possibly seen at 170 kpc north of the nu- 
cleus (Section l3.3t . Also the ghost bubbles seen to the South in 
the X-ray pressure maps jpabian et al1|200 6^ may be larger and so 
require more power than the present inner bubbles. If correct, the 
central source thus does vary by a factor of few on timescales of 
10^ yr. Deeper X-ray observations at larger radii than the current 
ones are needed to test this interesting possibility further. 

We conclude the discussion of sound waves by noting that typ- 
ical sound waves will be difficult to see in current data on other 
clusters. The X-ray surface brightness of the central 50 kpc of the 
Perseus cluster is more than twice that of another cluster. Large 
ripples, interpreted as weak shocks have been seen in the Virgo 
cluster ( Forma n et alj|2006l) and possible outer ripples ma y occur 
in A 2199 and and 2A 03354-096 jSanders & Fabiaj|2006i) . These 
may just be the peaks in the distribution of sound waves in those 
clusters from the more exceptional power episodes of their central 
engines. Like tree rings or ice cores for the study of geological his- 
tory, ripples in cluster cores offer the potential to track the past his- 
tory of an AGN for more than 10^ yr. Observations of a yet wider 
region are required to see where the ripples eventually die out. 



9.2 The nonthermal component 

We have reported on a hard X-ray component coincident with the 
radio minihalo (section|7](. It is difficult to support a thermal origin 
for this emission in terms of a ~ 16keV gas associated with the 
HVC. A simple energy argument against the thermal hypothesis is 
to note that the crossing time of the 200 kpc diameter inner core 
of the cluster at 3000 km s^' takes nearly 7 x 10^ yr. The radiative 
cooling time of 16 keV gas with a plausible density of 10^^ cm^-' 
is 10^'* yr, so the shocked gas would radiate only 0.7 per cent of its 
energy. With the luminosity of this component at 5 x lO^"' erg s^' 
we obtain a total injected energy of 2 x 10^' erg. This is the total 



kinetic energy of 2 x lO'' Mp, moving at 3000 km s^' . So an im- 
plausible large mass needs to be stripped from an apparently small 
galaxy in order to explain the hard component as due to shocking 
by the HVC. 

The hard component is therefore most readily interpreted as 
inverse Compton em ission from the minihalo. As discussed by 
ISanders etal1 ( l2005h this means that it must be well out of equipar- 
tition for the electrons and magnetic field with the electrons domi- 
nating the nontherm al pressure (as also deduced for the radio lobes, 
lFabianetal.li2002bl) . It is possible that the radio bubbles leak a 
smalJ3 fraction of their cosmic-ray electron content (and presum- 
ably protons or positrons for charge neutrality) into their surround- 
ings. These accumulate, losing their energy principally through in- 
verse Compton losses on the Cosmic Microwave Radiation. A half- 
power break in the power-law spectrum is then expected around 
10 keV in the hard X-ray flux if this process has, as expected, been 
continuing for more than a Gyr. This matches what is inferred of 
the observed spectrum. 

The presence of the nonthermal component increases the heat- 
ing within the inner 30 kpc if SP is raised proportionately, and 
also means that some di rect collisional heating of gas is possible. 
iRuszkowski et al.l l l2007h have since investigated this idea in more 
detail. 



9.3 The distribution of metals 

The enhanced metallicity in the core of the Perseus cluster shows 
a spatial distribution which extends to the S along the axis of old 
bubbles, as expected if the bubbles push and drag the gas around 
feoedieer et al."2007). What is particularly interesting here is evi- 
dence that the distribution is clumpy and also that the central drop 
is a real effect. 

The dumpiness and especially the sharp edges of at least one 
clump (Fig.|13ll likely require a magnetic field configured to prevent 
dispersal. The central metallicity drop is not easily explained as just 
due to some outer gas falling in to replace inner gas dragged out, 
since it has a significantly lower entropy. It could be accounted for 
if the metals are highly inhomogeneous on a small scale, with the 
hig her metgillicity gas whi ch has the shorter cooling time cooling 
out dMorris & Fabianl !2003). 

An interesting possibility raised by the high nonthermal pres- 
sure near the centre is that t he gas is made buoyant by the cosmic 
rays l lChandranll2004ll2005l) . In this picture the simple entropy in- 
ferred from just the gas temperature and density are insufficient to 
determine its behaviour. The gas cosmic-ray mixture becomes con- 
vectively unstable where the cosmic-ray pressure begins to drop 
steeply outward, leading to the gas overturning. The pressure drop 
is inferred to occur at about 30 kpc (Fig.l24t and is about the radius 
where the metallicity peaks (Fig. |9j. Such a turnover of gas may 
happen sporadically when the cosmic-ray density has built up for 
some time. 



9.4 The optical filaments 

Finally we briefly discuss t he Hg filaments. T hese radiate most 
of their energy at Lya (see iFabian et alj[l984l for an image) and 
are mostly composed of molecular gas at a few thousand down to 



^ Studies of bub bles in nearby clusters indica te they are not magnetic pres- 
sure dominated jPunn & Fab"i£u]|2004 1200^ . so this limits the number of 
particles which can escape. 
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50 K jHatch et alj2005l : |johnstone et al.l2007l : ISalome et alj2006h . 

We have shown that they are surrounded and mixed with soft X-ray 
emitting plasma at 0.5-1 keV temperature. In section|6]we consider 
one filament in detail. The Ha luminosit y of its p eak is about one 
per cent of the total measured by Heck man et al. 

I (1989), so scal- 
ing the X-ray inferred mass cooling rate there of 0.06 Mp yr^ ' we 
obtain a total mass cooling rate into filaments of 5 yr^' . 

This value assumes however that the X-ray emitting gas loses 
its energy solely by radiating X-rays. If in addition it loses energy 
by conduction, mixing or other means, with the cold gas in the 
filament and thereby powers the Lya emission then we need to 
scale the above rate by a factor of 20 (the expected Lya/Ha ra- 
tio for recombination) to obtain a total mass cooling rate of about 
100 Mp yr^' . This value agrees of course with that obtained by 
just taking the total Ha luminosity and assuming it is obtained 
from the thermal energy of the hot gas. It is about one third of 
the total inferred mass cooling rate obtained if there is no heating. 
We note that radiation at other wavelengths such as O VI emission 
teregman et alj200fl ) can increase this fraction. 

This implies that a high non-radiative mass cooling rate is pos- 
sible in the Perseus cluster. By extension, it suggest that this hap- 
pens in most cool-core clusters (which generally also have optical 
filaments; see Crawf ord et al. 1999). It can account f or part of the 
lack o f cool X-ray emitting gas in such cluster cores jpabian et al.l 
l2002ah . 

A cooling rate of lOOMQyr^' for 5 Gyr gives a total of 
5 X lO" Mq of cold gas. This is a bout 10 times more than is in- 
ferred from CO measurements ( S alome et al .l2006h . Star formation 
is another possible sink. NGC 1275 has long been known to have an 
A- star spectrum and e xcess blue light. UV imaging of NGC 1275 
by ISmith et al.l jl992h shows a lack of stars above 5M0, which 
means that any continuous star formation (which would have been 
only 20 My, yr^' ) must have ended about 50 Myr ago. Burst mod- 
els of star formation at rates up to many l OOMp yr~^ 100 M yr 
or so ago are consistent with the UV data jSmith et alJI 19921) . A 
picture in which gas accumulates through filaments and is then 
converted into stars sporadically on a 100 Myr timescale appears 
possible. 

If however much of the power in the filaments is due to sources 
other than the hot gas, cosmic rays or kinetic motions for example, 
then the above star formation estimate is an upper limit. 

9.5 Summary 

We have quantified the properties of the X-ray surface brightness 
ripples found in the core of the Perseus cluster and, assuming that 
they are due to sound waves, have determined the power propa- 
gated as sound waves. The power found in this way is sufficient to 
balance radiative cooling within the inner 70 kpc, provided that it 
is dissipated as heat over this lengthscale. This provides the quasi- 
isotropic, relatively gentle, heating mechanism required to prevent 
a full cooling flow developing. Ultimately the power is derived 
from the jets emitted by the central black hole. 

A hard X-ray component is confirmed and argued to be plau- 
sibly due to inverse Compton scattering by cosmic-ray electrons 
in the radio minihalo. The electrons may have leaked out of the 
radio lobes of 3C84 and now have a pressure comparable to the 
thermal pressure of the hot gas in the innermost 30 kpc. The lobes 
themselves appear to be devoid of any thermal gas unless its tem- 
perature is very high (50-100 keV). The cosmic-ray electrons are 
important in enhancing the heating and possibly also in changing 
the convective stability of the central 30 kpc. 



The X-ray data provide insight on the history of the past 
100 Myr of activity of the nucleus of NGC 1275. There are hints 
from the large ripples beyond 100 kpc, and from the the large 
Northern filament and the presence of many A stars, of a higher 
level of activity before that, a few 100 Myr ago. This can be tested 
by deep, high spatial resolution, observations of a wider region than 
covered out so far. 

We infer that a close balance between heating and cooling 
is established in the core of the Perseus cluster over the past few 
100 Myr. The average heating rate is, and has been, close to the ra- 
diative cooling rate, although there can be variations by a factor of a 
few on longer timescales. The primary energy source is the central 
black hole and jets; the energy is distributed by the sound waves 
generated by the inflation of the lobes. We suspect that this process 
is common to most cool core clusters and groups and is the mech- 
anism by which heating of the cool core occurs. It will however be 
difficult to verify observationally in those other objects since the 
X-ray surface brightness is so much lower. Only the extreme peaks 
in the distribution of activity will generally be detectable. 
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APPENDIX A: A DIRECT SPECTRAL DEPROJECTION 
METHOD 

Al Problems with existing spectral deprojection methods 

Results from fitting cluster data using the PROJCT model in XSPEC 
can be misleading. PROJCT is a model to fit spectra from several 
annuli simultaneously, to account for projection. There are one or 
more components per deprojected annulus, each with parameters 
(e.g. temperature, metallicity). The projected sum of the compo- 
nents along line of sights (with appropriate geometric factors) are 
fitted against each of the spectra simultaneously. Often the result- 
ing deprojected profiles (e.g. temperature) oscillate between values 
separated by several times the uncertainties on the values. This os- 
cillation can disappear if different sized annuli are used. Sometimes 
halving the annulus width can halve the oscillation period, indicat- 
ing they are not physical changes on the sky. 

The oscillation can be alleviated by fitting the shells sequen- 
tially from the outside, freezing the parameters of components in 
outer shells before fi tting spectra from shells inside them (see e.g. 
ISanders et al.ll2004h . This helps to solve the issue where poorly 
modelled spectra near the centre can affect the results in outer 
annuli (The standard way to fit the data is simultaneous. All the 
spectra are used to calculate each point, even though interior shells 
and not projected in front of outer shells.) The difficulty with this 
method is that uncertainties calculated on parameters to the model 
are underestimated. They do not include the uncertainties on outer 
shells. 

The outside-first fitting procedure does not fix every oscillat- 
ing profile. Numerical experiments, when clusters are simulated 
and fit with PROJCT (R. Johnstone, private communication), show 
that assuming the incorrect geometry when trying to account for 
projection does not produce oscillating profiles. Something which 
does create oscillating profiles are shells which contain more than 
one spectral component (e.g. several temperatures). It appears that 
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PROJCT tends to account for one of the components in one of the 
annulus fit results, and another in a different shell. By assuming that 
a spectral model is a good fit to the data, a very misleading result 
is produced. Other deprojection methods which assume a spectral 
model to do the deprojection will have similar issues. 

A2 Direct spectral deprojection 

We describe here a method to create 'dcprojcctcd spectra', which 
appears to alleviate some of the issues found using projct. It is 
a model independent approach, assuming only spherical geometry 
(at present). 

The routine takes spectra extracted from annuli in a sector, 
and their blank-sky background equivalents. From each of these 
foreground count rate spectra we subtract the respective blank sky 
background spectrum. Taking the outer spectrum, we assume that it 
was emitted from part of a spherical shell, and calculate the (count 
rate) spectrum per unit volume. This is then scaled by the volume 
projected onto the next innermost shell, and subtracted from the 
(count rate) spectrum from that annulus. After subtraction we cal- 
culate a spectrum per unit volume for the next innermost shell. We 
move inwards shell-by-shell, subtracting from each the calculated 
contributions from outer shells. This yields a set of deprojected 
spectra which are then directly fitted by spectral models. 

To calculate the tmcertainties in the count rates in the spec- 
tral channels in each spectrum we used a Monte Carlo technique. 
Firstly each of the input foreground and background spectra are 
binned using the same spectral binning, using a large number of 
coimts per spectra channel (we used 100 in this work) so that Gaus- 
sian errors can be assumed in each spectral bin. We repeat the de- 
projection process 5000 times, creating new input foreground and 
background spectra by simulating spectra drawn from Gaussian 
distributions based on the initial spectra and their uncertainties. The 
output spectra are the median output spectra from this process, and 
the 15.85 and 84.15 percentile spectra were used to calculate the 
1(7 errors on the count rates in each channel. 

This technique assumes that the response of the detector does 
not change significantly over the detector. This is the case for the 
ACIS-S3 detector on Chandra used here. It also assumes that the 
effective area does not change significantly, as we do not account 
for the variation of the ancillary response. This could be incorpo- 
rated into this method, but we have not implemented this yet. 
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Figure 16. Emission measure profiles across the filaments in the differ- 
ent temperature components. The top panel shows an unsharp-masked 0.5- 
7 keV X-ray image rotated so that the bins lie across it. The second panel 
shows a similar Ha image. The next panels show the 0.5, 1, 2 and 4 keV 
temperature component XSPEC normalisations, measured from the 1 arc- 
sec wide bins. The final panels show the best fitting absorbing Hydrogen 
column density and the metallicity of thermal components. 
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Figure 17. Results of fitting a cooHng flow model to spectra extracted across 
the filaments in Fig. 1151 The top panel shows an image of the region. In the 
second is plotted the mass deposition rates per bin obtained by fixing the 
lower temperature of the cooling flow models to certain values. The bottom 
panel shows the best fitting lower temperature if it is allowed to be free. 
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Figure 20. Flux per unit area on the sky for multicomponent mod- 
els plus powerlaws or plus a hot 16 keV component. Also pl otted are 
the b ackground-subh'acted results from an earlier analysis <Sanders et alj 
|2005|) . 
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Figure 18. Comparison o f Hg surface brightne ss flux profile across the 
filaments (using data from lConselice et"al]|200ll) versus unabsorbed 0.001 
to 30 keV surface brightness flux from the cooling flow model cooling down 
to 0.25 keV. Cooling to 0.0808 keV, the minimum available, increases the 
results by about 5 per cent. Both of these fluxes are in the same units. 
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Figure 21. Fluxes for the powerlaw or hot thermal components for various 
models. The units are logjg erg cm"^ s^ ' , in the 2-10 keV band. The panels 
show, from left to right, top to bottom, (1) fitting a thermal MEKAL plus a 
r = 2 powerlaw from 0.6 to 8 keV, (2) only fitting it between 2.3 and 8 keV, 
(3) allowing the T to vary from 1.4 to 2.4, (4) multicomponent thermal 
APEC model plus a F = 2 powerlaw, (5) using F = 1.8, (6) using F = 1.5, 
(7) allowing Fto vary between 1.4 and 2.4, and (8) multicomponent thermal 
APEC model plus 16 keV component. 
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Figure 22. Nw column density maps in units of 10^^ cm^^ generated by 
fitting different spectral models to regions containing a signal to noise ratio 
of 500. Also shown is an X-ray image of the same area. 'Single' shows the 
results of fitting a single PHABS absorbed MEKAL. 'Gamma=2, 1.8 and 1.5' 
shows the results fitting a multitemperature model plus a powerlaw of the 
photon index given. '16 keV shows the results using a multitemperature 
model including a 16 keV hot thermal component. 
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Figure 23. Spectrum showing the effect of the 16 keV component to spec- 
tral fit. The spectrum shown is a region around 1.8 ai'cmin north from the 
nucleus containing ~ 2.5 x 10^ counts. The solid line in the top panel is the 
model, after removing the 16 keV component. The bottom panel is the ratio 
of the data to the model not including the hard component. 
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Figure 24. Inferred average nonthermal particle pressure calculated from 
the r= 1.5 powerlaw plus multitemperature results in Fig. 1211 assuming 
inverse Compton emission. Also plotted is the average thermal gas electron 
pressure from Sanders et al. 1 2004). 
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Figure 25. Regions used for examining the spectra inside and outside of the 
bubbles. 
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Figure 26. Hot thermal component normalisation per unit area (defined by 
the BACKSCAL header keyword) as a function of temperature inside the 
bubbles compared to outside. The results of including a F = 2 powerlaw 
component are also shown. 
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Figure 27. 2a upper limits to the volume filling fraction of the bubbles. 
Models including a F = 2 powerlaw component are indicated. 
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